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The possibility of a friction term in the equation of motion for a scalar field is investigated in 
non-equilibrium field theory. The results obtained differ greatly from existing estimates based on 
linear response theory, and suggest that dissipation is not well represented by a term of the form 
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Most realizations of inflationary cosmology suppose 
that at some point in the universe's early history, its 
energy content became dominated by a classical scalar 
field (the inflaton 4>(t)), whose equation of motion has 
the form 



4> + 3(a/a)^ + r,(4>)^ + VM = 0, 



(1) 



where a(t) is the cosmic scale factor. The friction term 
rj{<j))(j) is important, because it provides a mechanism 
whereby the energy of the inflaton may be converted into 
the matter with which the present universe is filled. It 
has also been speculated |l| ^ that sufficiently large fric- 
tional effects might exercise a controlling influence on the 
inflationary process itself. To derive this term from fun- 
damental physics is, however, not straightforward. In 
fact, we argue in this letter that existing estimates of the 
coefficient r](tf>) are incorrect, and that dissipation prob- 
ably cannot be represented by a term r](j) at all. 

The most systematic derivations to be found in the 
literature rely on linear response theory || (see also 
@' H' B 01 )■ ^° ^ e concre te, suppose that cj)(t) is the 
expectation value of a quantum field <f>(t, x) with the self- 
interaction (A/4!)<I> 4 and write $(t, a?) = <p(t) + \(t,x). 
It is sufficient (see below) to deal with the Minkowski- 
space version of this theory, for which the exact equation 
of motion is 



= 0. 



(2) 



We hope to extract from the expectation value (x 2 (A x )) 
a contribution proportional to <f), and the ellipsis indi- 
cates terms that are irrelevant for this purpose. The 
particles associated with the quantum field x{t, x) have 



a time-dependent effective mass m (t) 



and, correspondingly, a time-dependent energy u>k(t) = 
\Jk 2 + m 2 (t), with 6j k = X(f>4> /2u> k . To use linear re- 
sponse theory, we assume that cj>(t) varies slowly on a 
time scale set by the relaxation times TjT 1 that charac- 
terize the evolution of the state of the x particles. If 
we wish to evaluate the equation of motion at time t, 
when tj> has the value tfio, say, then at times t' not too 
different from t we write <j>(t') = <fio + 6<p(t') and treat 
5<f>{t') ~ <j>(t)(t' — t) as a small correction. The Hamil- 
tonian for \ is H(t') rs Hq + 5H(t'), where Hq is the 



constant Hamiltonian obtained by setting cj>(t) — tfrg and 
SH(t') = {\/2)8<t>{t') J d 3 x [ X 2 + |x 3 ] • We obtain an ap- 
proximation to (x 2 (t,x)) of the Kubo type, namely 



<X*(t,x)) « ( X 2 (t,x)r + t / dt'([SH(t'),x 2 (t,x)]) 



(3) 

where the expectation values are taken in the equilibrium 
state associated with the Hamiltonian Hq, at a fixed tem- 
perature f3 . The second term, proportional to <fi, can 
be calculated by the methods of equilibrium perturba- 
tion theory j8j. To obtain a finite result, one must ap- 
proximately resum self-energy insertions, which have an 
imaginary part Tk ■ Assuming that Tk C^, the friction 
coefficient is 



\ 2 4>l f d 3 k 



(2tt) 3 



where n^. q (i) = [exp(/3uk(t)) — 1] 



(4) 

is the Bose-Einstein 



distribution function. 

Two heuristic arguments [|[ ^ indicate the origins of 
the terms in (Q). Morikawa and Sasaki Jl treat %(t, x) 
in the first instance as a free field. It can be expanded 
in terms of annihilation and creation operators ak (t) and 
at (t) , chosen to diagonalize the free- field Hamiltonian at 
time t, whose time dependence arises, in the usual Bo- 
goliubov formalism, from the time-dependent mass m(t) 
- an effect conventionally described as 'particle creation'. 
We define two functions rik(t) and Vk(t) through 



(4(*) a fe'W> 

(a k (t)a k > (t)) 



{2nf5{k - k')n k {t) 
(2n) 3 5(k + k r )v k (t). 



(5) 
(6) 



Of these, n k (t) (equal to n e k q (t) in a state of local equi- 
librium) counts the number of particles per unit volume 
with momentum k, these 'particles' being defined as the 
quanta annihilated by a k (t), while v k (t) measures the ex- 
tent to which the density matrix at time t is off-diagonal 
in the Fock basis associated with this definition of a par- 
ticle. The expectation value of interest is given by 



{X 2 (t,x))) 



d 3 k 



{2n) 3 2u k 



[l + 2n* + 2Rei*], (7) 



2 



while n k (t) and v k (t) obey the exact evolution equations 

(8) 



d t n k (t) = ^4irRe^(i) 



LO k (t) 



U! k (t) 



dtv k {t) = -2iu k (t)v k {t) + -^-[l + 2n k {t)]. (9) 

2u) k (t) 

It is argued in Ref. B that interactions introduce an 
imaginary part to the self-energy, in effect replacing tv k 
with io k — iT k . Then, with the further assumption that 
all the functions in (|9|) are approximately constant over 
a time interval of order T^ 1 , this equation can be inte- 
grated, yielding a contribution to v k {t) that reproduces 
the second term of (^), except that n k (t) need not have 
its equilibrium value. 

Hosoya and Sakagami @] identify the first term of (^) 
as arising from n k (t) in (M) by assuming that the evolu- 
tion of this function is described by the kinetic equation 



d t n k (t) = -2T k Mf)-< q (t)] 



(10) 



sometimes called the relaxation-time approximation to 



the Boltzmann equation. On writing n k (t) 
5n k lt), with dtn 



cq 



cq 



(*) 



-Pib k n k q (l + n k q ), and assuming 
again that everything changes sufficiently slowly on a 
time scale of IT , we obtain a solution for 5n k which, 
when inserted into (Q), reproduces the first term of (^J). 

The ad hoc features of these heuristic arguments 
(specifically, the fact that equations (||) and apply 
only to a free field theory, while the kinetic equation Jl(i| ) 
is essentially a guess) are largely avoided by the linear re- 
sponse calculations based on (||), developed by the same 
authors. However, the evaluation of the expectation val- 
ues in (|J) relics in an essential way on analytic proper- 
ties of thermal Green functions (arising from the so-called 
KMS condition) which apply only to a state of exact ther- 
mal equilibrium, and have no analogue for a system that 
departs, even to a small extent, from equilibrium. Since 
the equation of motion (Q) implies a non-equilibrium 
state in which <fi changes with time, it is important to 
know whether the result ([|) applies to the slow-evolution 
limit of a non-equilibrium state. If it does, it might also 
serve as a useful approximation when the evolution is not 
especially slow. (Even in equilibrium, the expression (Q) 
may be unreliable owing to large corrections from higher 
orders of perturbation theory [l^, 11 , but we do not 
address that issue here.) 

To investigate this question, we have derived improved 
versions of equations (||) - (|h]) which, though approxi- 
mate, take systematic account both of interactions and of 
non-equilibrium time evolution. They can be expressed 
as 

ed t n k = a k ~ T k (l + 2n k ) + e— Rev k (11) 

ed t v k = -2i{u k ~iT k )v k -a k + e^-{l + 2n k ), (12) 

2uj k 



where a k is a function whose meaning will be explained 
shortly, while e, which has the value 1, has been intro- 
duced as a formal parameter to generate a time-derivative 
expansion. 

Details of the somewhat lengthy derivation of these 
equations will be reported elsewhere, but the essential 
strategy is the following. We deal with Green functions 
p{t,x;t',x') = i(\x(t,x),x(t',x')}) and C{t, x; t', x') = 
h{{x(t,x),x(t' ,x')}}, where [. . .] denotes the commuta- 
tor and {. . .} the anticommutator. In particular, the ex- 
pectation value that we wish to evaluate is given by 



( X 2 (t,x)} = C(t,x;t,x). 



(13) 



These Green functions obey exact Dyson-Schwinger 
equations, which are conveniently expressed (see, e.g., 
Ref. [jl2| ) in terms of self-energies E p and T,q as 

[d* + u 2 k (t)}p(t,t';k) = - J dt"V p {t : t"-k)p(t",t':k) 

(14) 

[d 2 t +ul(t)]C(t,t'-k) = - ( dt'"E p (t,t";k)C(t",t';k) 



I dt'"E c (t,t";k)p(t",t';k) 
Jo 

(15) 



after a spatial Fourier transform. Our basic approxima- 
tion is to introduce local ansatze for the self-energies. 
First, the commutator function p(t, t'; k) reduces in ther- 
mal equilibrium to the spectral density, whose temporal 
Fourier transform can reasonably (though not exactly) be 
represented by a Breit-Wigner function, characterized by 
a thermal width T k . The corresponding non-equilibrium 
approximation is equivalent to the assumption 



E p (t, t'; k) n 2 [T k (t)T k (t')} 1/2 d t S(t - t' 



(16) 



Then, given that correlations decay roughly exponen- 
tially with time, on the time scale set by T k (t), it is 
reasonable to approximate £cs which is symmetric in its 
time arguments, as 



Eo(*, t'; k) « -2u> k (t)a k (t)5(t - t'). 



(17) 



The functions T k (t) and a k (t) remain to be determined. 
On substituting these ansatze into the Dyson-Schwinger 
equations, we obtain a pair of local differential equations 
for approximate Green functions, say p and C. These 
equations can be solved in terms of auxiliary functions 
n k {t) and v k (t) that satisfy the evolution equations ( p"l| ) 
and (|l2]). In particular, C(t,x;t,x) is then given pre- 
cisely by the expression (0). In the interacting theory, 
though, one cannot identify n k (t) and v k (t) unambigu- 
ously in terms of the expectation values (||) and (|^) . 

The approximate Green functions obtained in this way 
can be taken to serve as the lowest-order propagators in a 



3 



partially-resummed perturbation theory, similar to that 
described in Ref. The result of substituting the 

solutions of (|ll|) and (|l^) into the expression (Q) is the 
expectation value (x 2 (t,x)) evaluated at the lowest or- 
der of this perturbation theory. Here, we shall be satisfied 
with this lowest-order estimate, though the approxima- 
tion can in principle be pursued to arbitrarily high or- 
ders. Equally, this perturbation theory can be used to 
evaluate the self-energies in © and (0). With a suit- 
able prescription for extracting local approximations to 
these self-energies, we obtain concrete expressions for the 
functions a k (t) and Tk(t). At 2-loop order, we find 

c f A 
a k (t) = — / d 3 k 1 d 3 k 2 d 3 k 3 

71" J L0kUk 1 Uk 2 Uk 3 

x [(1 + n kl )(l + n k2 )n k3 + n kl n k2 (l + n k3 )] 

(18) 

r fc (t) = 4 / ^k^k^k 



(H), with which we plan to compare our own. With this 
approximation, we find that 



7T J U) k 0J kl U) k2 tjJ k3 

x [(1 + n fel )(l + n k2 )n k3 - n kl n k2 (l + n k3 )] 

(19) 

where c = A 2 7r 2 /64(27r) 5 and A is the product of delta 
functions that conserve energy and momentum in the 2- 
body collisions (k 1 ,k 2 ) # (k 3 ,k). 

With these expressions for a k and T k , the evolu- 
tion equation ( |Tl| ) for n k (t) is just what we might ex- 
pect: a genuine Boltzmann equation in place of the ap- 
proximate version hypothesized in (|l0|). The quantity 
S k = oi k — T k (l + 2n k ) is the scattering integral that 
accounts for the net rate of change of the number of par- 
ticles with momentum k due to collisions, while the re- 
maining term is the source term from (^) associated with 
particle creation. An heuristic expectation for the evo- 
lution of v k (t) is not easy to formulate. We can observe 
in (|l2]), though, that u> k has indeed been replaced with 
uj k — iT k in the first term of (JsJ) , though not elsewhere. 

We investigate the solution of (11 ) and ( |l2] ) in the limit 
of slow time evolution by expanding in powers of the 
formal parameter e, so as to generate a time-derivative 
expansion. For n k (t), the leading term is just the Bosc- 
Einstein distribution n c k q (t), which makes S k vanish. The 
inverse temperature [3 is, of course, a free parameter spec- 
ifying a particular state. For v k (t), the leading term is 



ia k q {t) 



2[u k {t)- l T cq (t)] 2[uk{t)-iTl\t)] 



(20) 



where ct c k q (t) and T° k q (t) are the above expressions with 
n k = n° k q (t). The next-to-leading terms Sn k (t) and 5v k (t) 
are proportional to <fr, and it is from these that we can 
obtain rj(<p) by substitution in (^). Here, we use the aux- 
iliary approximation that T° k q <C u>k- This slightly simpli- 
fies the numerical calculations to which we must shortly 
resort, but we have adopted it primarily for consistency 
with the strategy that leads to the linear-response result 



Re5v k = - — 2 
2uji 



^rl q (l + 2nt q )-u k d t Imvl q 



while Sn k is the solution of the integral equation 



dk'K(k, k')5n k , = d t n q - —Rev cq , 



(21) 



(22) 



with K(k,k') = SS k /Sn k >\ nk=n ^. 

At this point, we observe that the linear-response re- 
sult can be recovered by means of two further approx- 
imations. The first is to set K{k, k') — —2T^5(k — 
k'), which means replacing a k and T k in (O) with 
their equilibrium values. The second is to take the 
time dependence of T k (t) to be of the form T k (t) « 
const. /u>k(t). This amounts to replacing the ansatz ( |l6| ) 
with T k (t) iY] c p q (uj k ,k)/2uj k (t), and taking the equilib- 
rium self-energy to be constant. 

To avoid these extra approximations, we must, in par- 
ticular, solve the integral equation (|2|) using the full ker- 
nel K(k,k'). This equation is in fact not self-consistent. 
The reason is that the scattering integral S k conserves 
both energy and (at 2-loop order where only elastic col- 
lisions are included) particle number. This implies two 
sum rules 



(23) 



/ dkk 2 K(k,k')= dkk 2 uj k K(k, k') = 0, 
Jo Jo 

which are not respected by the right-hand side of (|22 
In principle, the equation has no solution. This means 
that the time-derivative expansion does not exist, and 
that dissipation cannot properly be represented by a term 
proportional to <f> in the equation of motion. Neverthe- 
less, we have obtained numerical solutions to ( f22| ) that 
are rather well defined. This is possible because the 
right-hand side becomes very small as k increases be- 
yond quite modest values and the same turns out to be 
true of the solution Sn k . In effect, the solution involves 
only a truncated kernel, which is not constrained by the 
sum rules. The existence of this approximate solution al- 
lows us to construct an effective friction coefficient, and 
thus to quantify the effect of the extra approximations 
that lead to (f|). The effect is rather large, as indicated 
in figures [l] and 0. Figure [l] shows the dimensionless 
quantity a = 77(0)rn/(1287r0 2 ), calculated from the lin- 
ear response formula (Q) as a function of the coupling 
strength c defined above, for selected values of the di- 
mensionless parameter /3m. The haphazard appearance 
of the four curves, which have essentially the same shape 
when viewed on appropriate scales, is explained by a com- 
petition between the two terms in (Q). Since T k is propor- 
tional to c, weak coupling favours the first term, which is 
a decreasing function of temperature, while strong cou- 
pling favours the second term, which increases with tem- 
perature. Figure ^| shows a as calculated numerically 



4 




FIG. 1: The friction coefficient a as calculated from the linear 
response formula (^). 



pm = 5.0 
pm = 1.0 
pm = 0.5 
pm = 0.1 




ond term in (g) where, for the mode of momentum k, the 
time integral amounts to an average over a time interval 
of order FT . During this time interval, the actual state 
of the system changes by an amount of order T^cj), so 
the putative equilibrium state represented by (. . .) oq is 
uncertain by an amount of this order. Thus, the first 
term of (|J) is ill-defined by an amount proportional to 4> 
that could well be comparable with, or greater than, the 
quantity we extract from the second. 

Two final remarks are in order. First, the field the- 
ory in an expanding universe is equivalent, in confor- 
mal time, to a Minkowski-space theory with a mass m$ 
that depends on a(t). Applied to this theory, the time- 
derivative expansion generates additional terms, propor- 
tional to a/a, which may be important, but do not affect 
our conclusions about the term proportional to <fr. Sec- 
ond, the results presented here may be specific to the A<& 4 
model. The dynamics of fermions to which $ might cou- 
ple can be re pre sented fll4| by evolution equations similar 
to (0) and @, but K(k,k') and T k are different, and 
the solution may be very sensitive to this difference. 

It is a pleasure to acknowledge discussions with Arjun 
Berera and Rudnei Ramos, which prompted the study 
reported here, and to thank Jim Morgan for valuable 
comments on an earlier version of this paper. 



FIG. 2: The friction coefficient a as calculated numerically 
from (El) and (Eh. 



from equations (^l]) and (p2[). Not only does the mag- 
nitude differ by factors of several hundred, but a turns 
out to be negative at weak coupling. (Strictly, it is only 
at weak coupling that the assumption Tk <C Wfc holds.) 
Of course, a negative friction coefficient does not make 
good physical sense: it clearly reflects the breakdown of 
the time-derivative expansion noted above. 

We estimate the accuracy of the numerical integrals 
used to obtain K(k, k') and to evaluate (0) as about 1%, 
and have verified the sum rules (^3|) at this level. The so- 
lution of (|2^) for Snk is such that the left and right-hand 
sides of this equation differ by no more than one part 
in 10 7 (i.e. the solution is valid to full single-precision 
accuracy) and we have verified that the whole procedure 
converges to the results shown as the momentum cutoff 
used in (Q) and (^||) is increased. 

These results strongly indicate that linear response 
theory gives the wrong answer to this problem, and the 
reason may be this. The formula (||) arises from the sec- 
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